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Abstract 

Interactions between non-BPS non-Abelian vortices are studied in non-Abelian U(l) x 
SU(N) extensions of the Abelian-Higgs model in four dimensions. The distinctive feature 
of a non-Abelian vortex is the presence of an internal CP^" 1 space of orientational degrees 
of freedom. For fine-tuned values of the couplings, the vortices are BPS and there is no 
net force between two static parallel vortices at arbitrary distance. On the other hand, for 
generic values of the couplings the interactions between two vortices depend non-trivially on 
their relative internal orientations. We discuss the problem both with a numerical approach 
(valid for small deviations from the BPS limit) and in a semi-analytical way (valid at large 
vortex separations). The interactions can be classified with respect to their asymptotic 
property at large vortex separation. In a simpler fine-tuned model, we find two regimes 
which are quite similar to the usual type I /II Abelian superconductors. In the generic 
model we find other two new regimes: type I* /II*. Unlike the type I (type II) case, where 
the interaction is always attractive (repulsive), the type I* and II* have both attractive and 
repulsive interactions depending on the relative orientation. We have found a rich variety 
of interactions at small vortex separations. For some values of the couplings, a bound state 
of two static vortices at a non-zero distance exists. 
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1 Introduction 



Nobody doubts the importance of topological solitons in various areas of modern physics (see 
[1] for a general review). They are closely related to the phenomena of spontaneous symmetry 
breaking which e.g. occur as a phase transition from the high temperature phase of the early 
universe to the present cold universe. In particular, vortex strings are believed to play important 
roles in the confinement of quarks in QCD, and they could be relevant to the study of cosmic 
string effects in the early universe. Historically, the vortex string, as a topological soliton, was 
found in the Abelian-Higgs model by Abrikosov and Nielsen-Olesen [2} [3]. 

Recently, a new type of vortex was found in U(N) non-Abelian gauge theories coupled with 

Nf = N Higgs fields in the fundamental representation [HE]. This object is called non-Abelian 

vortex. A typical feature of this vortex is that it possesses internal degrees of freedom, which 

arise when the vortex breaks an exact flavor symmetry of the vacuum. These degrees of freedom 

are related to the orientations of the non-Abelian flux inside the vortex core. The vacuum of 

the theory leaves the colour-flavour locked SU(N)c+f symmetry; on the other hand, the vortex 

soliton spontaneously breaks this SU (iV) c+F symmetry to SU (N — 1)c+f x U (1)c+f; these broken 

symmetries give rise to the moduli space of an elementary vortex 

cpN -i = SU(N) C+F 

SU(N-l) c+F xU(l) c+F 

The classical moduli coordinate can be promoted to a field living in the vortex worldvolume; in 

this way vortex solitons in a 3 + 1 dimensional theory are directly connected with a CP N ~ l sigma 

model in 1 + 1 dimension, which describes the macroscopic physics of the flux tube. Several 

groups have intensively investigated these objects in relations to various aspects of physics; a 

partial list includes confined monopoles [6], quantum aspects [HE], [9], higher winding numbers 

[TUl [TT| [12], relation to D-branes in string theory [U [13], dualities [TU [151 US], cosmic strings 

[TTl UH] , semilocal extensions [T§1 [20], SO(N) generalization [21], high temperature QCD [2"2] . 

global vortices [23], gravity [21], composite states of various BPS solitons O [25] and statistical 

mechanics [26]. Readers can find good reviews in [27[ [28]. 

Most the works on the non-Abelian vortex so far were focused on the BPS limit [29] (a 
single non-Abelian non-BPS vortex configuration is discussed in [30l 1311 [T^ 122]). No forces arise 
among BPS vortices, because there is a nice balance between the repulsive forces mediated by the 
vector particles and attractive forces mediated by the scalar particles. In this particular limit, 
the solutions to the equations of motion develop a full moduli space of solutions [32]. However, 
once the balance between the attractive force and the repulsive force is lost, the moduli space 
disappears. Alternatively we can think that an effective potential is generated on this moduli 
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space. It is well known that ANO vortices in the type I system feel an attractive force while 
those in the type II model feel a repulsive force [291 E31 Ell ESI E3 EH] • In condensed matter 
physics, it is also known that type II vortices form the so-called Abrikosov lattice, [21 EH] due 
to the repulsive force between them. Furthermore, lattice simulations give some evidence of the 
presence of a (marginal) type II superconductivity in QCD [40J. 

We are interested in studying interactions between non-Abelian vortices which are non-BPS. 
In non-supersymmetric theories, BPS configurations are obtained with fine-tuned values of the 
couplings. If supersymmetry exists in the real world, it is surely broken at a low energy scale; 
therefore non-BPS vortices are more natural than BPS ones. Also, we encounter such non-BPS 
configurations in supersymmetric theories [51 [121 CE] when we consider a hierarchical symmetry 
breaking closely related to a dual picture of color confinement of truly non-Abelian kind. Specif- 
ically we are interested in the interactions between vortices with different internal orientations, 
which is the distinct feature from the ANO case. In a previous paper [Jl], we have discussed 
these aspects in an M = 2 theory with an adjoint mass term which breaks the extended super- 
symmetry, and we have found a natural non-Abelian generalization of type I superconductors. 
Even if the force between two non-Abelian vortices is not always attractive, we have found a 
close resemblance with type I Abelian vortices: the lightest field of the theory is a scalar field. 
So if we put two vortices at large distance, the prevailing part of the interaction is mediated 
by the scalar particles and not by vector particles. Moreover, if the two vortices have the same 
orientation in the internal moduli space, the force is always attractive. 

In this paper we study the same problem in another theoretical setting: an extension of the 
Abelian-Higgs model with arbitrary scalar couplings which is generically incompatible with the 
BPS limit. The simplest extension in this direction is a theory in which there are just two mass 
scales; the mass of the vector bosons and the mass of the scalars. There is one parameter A which 
controls the ratio of the two mass scales. We find that A < 1 leads to an attractive force as a 
usual Abelian type I, while for A > 1 a repulsive force works, similarly to the usual Abelian type 
II. There is no force between vortices with opposite CP N_1 orientation. For A < 1 (type I) this 
configuration is unstable and the true minimum of the potential corresponds to two coincident 
vortices with the same orientation. For A > 1 (type II) this configuration is stable; in other words 
a part of the moduli space corresponding to the relative distance between vortices with opposite 
orientations survives the non-BPS perturbation. 

However, in more general theories where the masses of the Abelian and non-Abelian degrees 
of freedom are different, we find a more complicated picture. There are four mass scales, the 
masses of the U(l) and of the SU(N) vector bosons, the masses of the scalars in the adjoint 
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representation and the singlet of SU(N)q+f- At large distance, the interaction between two 
vortices is dominated by the particle with the lightest mass. So if we keep the four masses as 
generic parameters, at large vortex separation we find four different regimes that we call Type 
I, Type II, Type I* and Type II*. In the last two categories repulsive and attractive interactions 
depend on the relative orientation. We study also numerically the interactions among two vortices 
at any separation with arbitrary orientations, and find that short distance forces also have rich 
qualitative features depending both on the relative orientations and the relative distance. 

The paper is organized as follows. In Sect. 2 we describe the theoretical set-up. In Sect. 3 we 
write the equations for the vortex and we quickly review the moduli space of the two vortices in 
the BPS limit. In Sect. 4 vortices in a fine-tuned setup are studied; the effective vortex potential 
in the case of small deviations from the BPS limit is found numerically. In Sect. 5 a more general 
set-up with four independent parameters is discussed in the same way. In Sect. 6 the effective 
potential at large vortex separation is found using a semi-analytical approach. Sect. 7 contains 
the conclusions. In the Appendix we provide the link between the formalism of this paper and 
that of the companion paper [4"T] . 



2 Theoretical set-up 
2.1 A fine-tuned model 

Our natural starting point is the following non-Abelian, U(N), extension of the Abelian-Higgs 
model in four dimensions: 



£ = Tr 



1 \2 2 

—^F^ + V,H(V»Hy - -A (v 2 l N - HW) 2 



(2-1) 



Here, for simplicity we take the same gauge coupling g for both the U(l) and SU(N) groups, 
while X 2 g 2 /4 is a scalar coupling and v (> 0) determines the Higgs VEV. In this simple model 
we have only three couplings (g, X,v). The N by N matrix field H embodies N Higgs fields in 
the fundamental representation of U(N). There is also an SU(N) flavor symmetry which acts 
on H from the right hand side. The vacuum of the model is given by: 

ifift = v 2 l N . (2.2) 

The vacuum breaks completely the U (N) gauge symmetry, although a global color-flavor locking 
symmetry SU(N)q+f is preserved 

H -> U G HU F , U G = Ul U G e SU(N) G , U f g SU(N) f . (2.3) 
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The trace part TtH is a singlet under the color-flavor group and the traceless parts are in the 
adjoint representation. We have two mass scales, one for the vector bosons and the other for the 
scalar bosons. The U(l) and the SU(N) gauge vector bosons have both the same mass 

M u(l) = M SU {N)= gv. (2.4) 

The masses of the scalars are given by the eigenvalues of the mass matrix. We start with 2N 2 
real scalar fields in H: N 2 of them are eaten by the gauge bosons (the Higgs mechanism) and 
the other N 2 (one singlet and the rest adjoint) have same masses 

M s = M ad = Xgv. (2.5) 

When we choose the critical coupling A = 1 (BPS), the mass of these scalars is the same as the 
mass of the gauge bosons and the Lagrangian allows an M = 2 supersymmetric extension. The 
BPS vortices saturating the BPS energy bound admit infinitely degenerate set of solutions. 



2.2 Models with general couplings 

A straightforward generalization of the fine-tuned model (12. ip is to consider different gauge 
couplings, e for the U(l) part and g for the SU(N) part, and a slightly more general scalar 
potential 



C = Tr 



- V, (2.6) 



2g 2 ^ 2e 2J ^ J 

where we have defined = F A T A and f„ u = F°„T° with Tr(T A T B ) = 5 AB /2 and 

The scalar potential is: 

\ 2 2 N 2 — 1 , 2 2 

A=X 

= ^ Tt x2 + ^ Ti (x T°-v 2 1 n )\ (2.7) 

where 

N 2 -l N 2 -l N 2 -l 

X = HH^ = X°T° + xA T A i X =Y1 XATA = 2 Y1 { Hl] T AH i) tA - ( 2 - 8 ) 

A=l A=l A=l 

The Lagrangian has the same symmetries as the previous fine-tuned model (12. ip . The potential 
in Eq. (12.71) is the most general gauge invariant quartic potential which can be built with the 
matter content of the theory. The U(l) and the SU(N) vector bosons have different masses 

M u{ i) = ev, M SU (N) = gv. (2.9) 



1 In the case g = e it is more compact to use F. IV = + /„„. 
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Moreover, the singlet part of H has a mass M s different from that of the adjoint part M a d 



M s = X e ev, M ad = X g gv. (2.10) 

When we take equal couplings, g = e and A = A e = A 9 , the scalar potential reduces to the simple 
potential V g=e = ^j^Tr (X — v 2 1tv) 2 . For the critical values A e = X g = 1, the Lagrangian allows 
an M = 2 supersymmetric extension and then the model admits BPS vortices which saturate the 
BPS energy bound. 



3 Non-Abelian vortices in the fine-tuned model 

3.1 Vortex equations 

We study the fine-tuned model (12.11) through out this section. For convenience, let us make the 
following rescaling of fields and coordinates: 

gv 

The Lagrangian then in Eq. (12.11) takes the form 



H^vH, W,^gvW„ x^-^. (3.1) 



C 

C = = Tr 



(3.2) 



and the masses of vector and scalar bosons are rescaled to 

M u{1) = Msu(N) = 1, M s = M ad = A. (3.3) 

As explained in the introduction, the model with A < 1 (A > 1) in the Abelian case (N = 1) is 
called type I (type II) and the forces between vortices are attractive (repulsive). At the critical 
coupling A = 1, there are no forces between vortices, so that multiple vortices stably coexist. 

In order to construct non-BPS non-Abelian vortex solutions, we have to solve the following 
2nd order differential equations, derived from the Lagrangian (12.11) . 

V^ v - % - [H(V V H) ] - (V V H)H^] = 0, (3.4) 
A 2 

V^H + — (1 - HH ] ) H = 0. (3.5) 

From now on, we restrict ourselves to static configurations depending only on the coordinates 
x 1 ^ 2 . Here we introduce a complex notation 

, n n <9i - id 2 W x - iW 2 £>i - iV 2 , N 

z = x i +ix ^ d= , W = ^— -, V=^— — 2 - = d + iW. (3.6) 

Zi Zi Zi 
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The equation of motions are of course not gauge invariant but covariant. It might be better to 
study gauge invariant quantities instead of dealing with the original fields H and W^. For this 
purpose we rewrite our fields as follows 

W(z, z) = -iS- x {z, z)dS(z, z), H(z, z) = S~\z, z)H(z, z), (3.7) 

where S takes values in GL(N, C) and it is in the fundamental representation of U(N) while 
the gauge singlet H is an N x N complex matrix. There is an equivalence relation (S, H) ~ 
(V(z)S, V (z)H), where V(z) is a holomorphic GL(N,C) matrix with respect to z, because dif- 
ferent elements in the same equivalence class give us the same physical fields as in Eq. (13. 7p . The 
gauge group U (N) and the flavor symmetry act as follows 

S(z, z) -> U G S(z, z), H (z) -> H (z)U F , U G G U(N) G , U F G SU(N) F . (3.8) 

In order to write down the equations of motion (13.41) and (13.51) in a gauge invariant fashion, 
we introduce a gauge invariant quantity 

n(z,z) = S(z, z)S(z, z) ] . (3.9) 

With respect to the gauge invariant objects Q and H, the equations (13.41) and (I3.5P are written 
in the following form 

Ad 2 (fMT 1 ) - Hd ^Q- 1 ^} + dHH^Q- 1 = 0, (3.10) 

no (n- l BH^j + 8 (no + -r (p - hh^ n^H = o. (3.11) 

Notice that Eq. (13 . 1 0[) is a 3rd order differential equation. This is the price we have to pay 
in order to write down the equations of motion in terms of gauge invariant quantities. These 
equations must be solved with the following boundary conditions for k vortices: 

detH^z k , Q — > HH\ as z -> oo. (3.12) 

The field strength is given by 

F 12 = 2S~ 1 B (ttdtt- 1 ) S. (3.13) 

Notice that Eq. (13. lOj) is invariant under the SU(N) flavor symmetry while Eq. (13.111) is covariant. 
This leads to Nambu-Goldstone zero modes for vortex solutions. 
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3.2 BPS Limit 



To see the relation with the BPS equations, let us take a holomorphic function H with respect 
to z as 

H = H (z). (3.14) 
Then the equations (13.101) and (13. lip reduce to 

3 

d (ndQ- 1 ) + ^ (l - P^IT 1 ) = 0. (3.16) 



Ad (fMT 1 ) - HqHIq- 1 ] = 0, (3.15) 
A 2 



These two equations are consistent only in the BPS limit A = 1. The equation (13 . 1 6[) is the 
master equation for the BPS non-Abelian vortex and the holomorphic matrix Hq(z) is called 
the moduli matrix (TUJ EE]- For any given moduli matrix H (z), given the corresponding solu- 
tion to the master equation, the physical fields and H are obtained via Eq. (13. 7p . All the 
complex parameters contained in the moduli matrix are moduli of the BPS vortices. For exam- 
ple, the position of the vortices can be read from the moduli matrix as zeros of its determinant 
det Ho(zi) = 0. Furthermore, the number of vortices (the units of magnetic flux of the configura- 
tion) corresponds to the degree of det H (z) as a polynomial with respect to z. The classification 
of the moduli matrix for the BPS vortices is given in Ref. [10], [28] . 

From now on, we consider a U(2) gauge theory (N = 2), which is the minimal model for 
non-Abelian vortices. The minimal winding BPS vortex is described by two moduli matrices 

\ -V 1 J \0 z-z J 

The complex parameter z corresponds to the position of the vortex while the other parameters, 
b and b', parameterize the internal orientation. This modulus gives rise to an internal moduli 
space CP 1 p^Ol [28] . In fact, they are inhomogeneous coordinates for CP 1 and are related by the 
transition function b = 1/b'. This orientational modulus can easily be understood from a simple 
argument related to the symmetry of the theory as we will see below. 

A rigorous way to define the orientation of the non-Abelian vortex is to identify it with the 
null eigenvector of Hq(z) at the vortex position z = zq. For the moduli matrix in Eq. (I3.17P the 
orientational vectors are 





(3.18) 



Here "~" stands for an identification up to complex non zero factors: ~ \<f), A G C*. One can 
easily find a direct relation between the parameter b (b') and the broken SU(2)c + p. Since H 
transforms as Ho — > HqUy under the color-flavor group, the orientational vector transforms as 
— > £/p0. If we start with = (1,0) T , we can recover = (l,b') T by use of the color-flavor 
rotation as 

- CH;r) (in;,-) <•■•> 

with |«| 2 + \(3\ 2 = 1. Thus we identify b' and (3* /a* . In what follows, we will call two non-Abelian 
vortices with equal orientational vectors parallel, while when they have orthogonal orientational 
vectors we will call them anti-parallel. The reader must keep in mind that vortices are always 
parallel in real space. Throughout this paper, we use the words parallel and anti-parallel only 
referring to the internal orientational vectors. 

Generic configurations of two vortices at arbitrary positions and with arbitrary orientations 
are described by the moduli matrices [TUl [28] : 

The superscripts label patches covering the moduli space. One more patch similar to (2, 0) is 
needed to complete the full moduli space [101 [28]. The positions of the vortices are the roots 
of zf — (0 + 0)zj + 00 — rjfj — zf — a'zi — f3' — 0. By using translational symmetry we can set 
z i + z 2 — (0 + = a' = 0) without loss of generality. The orientation vectors are 0i = 
(?7, z\ — 0) and 2 = (77, ^2 ~ 0) for the (1, 1) patch, while they are 4> 1 = (1, a' z\ + b 1 ) 
and 2 < ' 2 ' ^ = (I? a/ -22 + fo r the (2,0) patch. Overall complex factor does not have physical 
meaning, so that each vector takes value on CP 1 . We can describe anti-parallel vortices only in 
the (1,1) patch when t] = fj, because of 0{02 = 0. On the other hand, we can describe parallel 
vortices only in the (2, 0) patch when a' Z\ + b' = a' Z2 + b' . 

For convenience, let us take a special subspace where fj = in the (1, 1) patch: 

H^J d Z ~ Z ° ~ V ) , z = z 1 = -z 2 = <l> = -l (3.21) 
V z + z J 

One can always recover generic points in Eq. (13.201) using flavor rotations. The parameters 
{zo,ri) and (/3,a',b') are related by the following relations: f3' = Zq, a' = I/77 and b' = —Zo/r). 
The orientational vectors are then of the form 



(1,1) I _ 1 r (1,1)1 




'1 \z=z — 1 - I ' ^2 




(3.22) 



When the vortices are coincident, however, the rank of the moduli matrix at the vortex 
position reduces. In this case we can no longer define two independent orientations for each 
vortex but we can only define an overall orientation. In fact, one can see that when zq = 0, the 
two orientations in Eq. H 3 . 2 2 f) are both equal to </> = (1, 0) T . We cannot really give to the 
parameter rj an exact physical meaning of a relative orientation between two coincident vortices. 
It is better, in this case, to consider this parameter merely as an internal degree of freedom of 
the composite vortex. When we take correctly into account both the parameter rj and the global 
flavour rotations that we previously factorized out, we recover the full moduli space for coincident 
vortices, which is WCPf 21l) ~ CP 2 /Z 2 [HI HH [12]. 

The definitions of the position and orientation of a single vortex can be rigorously extended 
to the non-BPS case by merely replacing H (z) with H (z,z). For configurations with several 
vortices, all the flat directions that are not related to Goldstone modes or translational symmetries 
will disappear. It is possible to use these definitions as constraints on the Hq(z, z) matrix to fix 
positions and orientations. Then our formalism allows us to study the static interactions of 
non-BPS configurations. 



4 Vortex interaction in the fine-tuned model 

We now concentrate on the fine-tuned model (13.21) . We will first calculate the masses of a special 
class of non-BPS coincident vortices. Then we will derive an effective potential for coincident 
almost BPS vortices but with generic value of the internal modulus parameter. Finally, we will 
compute an effective potential for two almost BPS vortices at any distance and with any relative 
orientations. 

4.1 (&4, k,2) coincident vortices 

The minimal winding solution in the non-Abelian gauge theory is a mere embedding of the ANO 
solution into the non-Abelian theory. This is obvious also from the moduli matrix view point. 
In fact, in the non-Abelian moduli matrix (13.171) we can put b (or b') to zero with a global 
flavour rotation. Henceforth we can recognize the moduli matrix for the single ANO vortex: 
Hq NO (z) = z — z as the only non-trivial element of the moduli matrix. 

This kind of embedding is also useful to investigate a simple non-BPS configuration. Let us 
start with the moduli matrix for a configuration of k coincident vortices. Since we have an axial 
symmetry around the k coincident vortices, we can make the following reasonable ansatz for Q 
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and H 



1 o 

w(r) 



#(o,D 



1 

f{r)z h 



(4.1) 



Note that f(r) = 1 means H^ ^ = Hq\z) which is nothing but the BPS solution. We will call 
the multiple vortex which is generated by the ansatz in Eq. ( 14.11) "(0, fc)-vortex". In terms of the 
two fields w(r) = e y< - r ) and f{r) the equations (13.101) and (13.111) reduce to the following form 



Y" + ly/ _ 2 ( f 1 + 
r /\ 



(/" + 1 ~ ± ^f - Y'fj - A 2 (1 - r 2fc / 2 e" y ) = 0, 



Y'" + -Y" - ^Y + e~ Y r 2k - l f 2 (2k - rY) = 0. 



The boundary conditions are 



(4.2) 
(4.3) 



F->2Hogr, Y'^2k/r, / -> (r -> oo) 
y -> 0, /' -> (r -> 0). 



(4.4) 
(4.5) 



Although it is quite impossible to solve these differential equation analytically, we can solve them 
numerically. The results for the single vortex (k = 1) are shown in Fig. [TJ where we used several 
different values of A . 
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Fig. 1: Numerical plots of Y (left) and / (right) for the single vortex: black for A = 1, red for A = 1.7, 
blue for A = 0.5. The broken line is 21ogr. 



When k > 2 it is possible that the ansatz (14.11) does not give the true solution (minimum of 
the energy) of the equations of motion (13.101) and (13. lip . This is because there could be repulsive 
forces between the vortices. With ansatz (14. ip we fix the positions of all the vortices at the 
origin by hand. The reduced equations (14 .2p and (14.31) are nevertheless still useful to investigate 
the interactions between two vortices. The results are listed in Table [TJ For A = 1, the masses 
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A 


k = 1 


k — 2 


k = 3 


0.6 


0.81305 


1.52625 


2.21205 


0.7 


0.86440 


1.65337 


2.42101 


0.8 


0.91231 


1.77407 


2.62115 


0.9 


0.95737 


1.88936 


2.81382 


1 


1.00000 


2.00000 


3.00000 


1.1 


1.04053 


2.10655 


3.18045 


1.2 


1.07922 


2.20944 


3.35575 


1.3 


1.11626 


2.30905 


3.52639 


1.4 


1.15182 


2.40566 


3.69276 



Table 1: Numerical value for the masses of coincident vortices. 



are identical to integer values, up to 10~ 5 order, which are nothing but the winding number 
of the vortices. Furthermore, our numerical results for generic A are in perfect agreement with 
the numerical value for ANO vortices obtained about 30 years ago by Jacobs and Rebbi [33] . 
As mentioned, this happens because the (0, /c)-vortex is obtained by embedding of the k ANO 
vortices. 

There is another type of composite configuration which can easily be analyzed numerically. 
These configurations are generated by the following ansatz for Q and H 

n uW-iM ° \ fl™ = f AM** ° V (4 . 6) 

V w 2 (r) J \ f 2 (r)z k * ) 

This ansatz corresponds to a configuration with k\ composite vortices which wind in the first 
diagonal U(l) subgroup of U(2) and with k 2 coincident vortices that wind the second diagonal 
£7(1) subgroup. The two sets of vortices can be considered each as embedded ANO vortices for 
the two decoupled Abelian subgroups. We refer to these decoupled non-Abelian vortices as a 
"(fci, A^-vortex". The mass of a (k\, A^-vortex is thus the sum of the mass of the (ki, 0)-vortex 
and that of the (0, A^-vortex. For example, the mass of (1, l)-vortex is double of the mass of the 
(0, 1)- vortex listed in the first column of Table [TJ As in the previous case, we get the minima of 
the energy under the constraint that the vortices are coincident. 

Because our fine-tuned non-Abelian model is a simple extension of the Abelian-Higgs model, 
we expect similar behavior for the interactions. Actually we have only one parameter A. Thus 
we will call the non-Abelian vortices for A < 1 type I, while they will be called type II for 
A > 1. From Fig. [21 in which is summarized the relevant data of Table [TJ we can argue which 
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(0, 2)-vortex 



, , . (1, l)-vortex 

A = 1: 2x single vortex 



. " (1, l)-vortex 



(0, 2)-vortex 

(a) type I (A < 1) (b) type II (A > 1) 

Fig. 2: Spectrum of the (0,2) and (1, 1) coincident vortices. 



kind of interaction appears between two non-Abelian vortices. In the type I case, the (0, 2)- 
vortex is energetically preferred to the (1, 1)- vortex, while in type II case the (1, l)-vortex is 
preferred. If the two vortices are separated sufficiently, we can ignore any interaction between 
them. Regardless of their orientations, the mass of two well separated vortices is twice that of 
the single vortex. This mass is equal to the mass of the (1, l)-vortex. Furthermore, it seems that 
the two separated anti-parallel vortices do not interact, and the energy does not depend on the 
relative distance. For the type I case, Fig. [2] suggests that the configuration with (2, 0)-vortices 
is the true minimum of the system. This means that attractive forces appears between vortices 
with different orientations. An attractive force also works in the internal space, which aligns 
the orientations. In the type II case, it seems that we do not have an isolated minimum of the 
energy. In fact, all the anti-parallel configurations with arbitrary distance have the same value 
of the energy. In the following sections we will confirm the picture we have outlined here. 

4.2 Effective potential for coincident vortices 

The dynamics of BPS solitons can be investigated by the so-called moduli approximation [52] . 
The effective action is a massless non-linear sigma model whose target space is the moduli 
space. The sigma model is obtained by plugging a BPS solution into the original Lagrangian 
and promoting the moduli parameters to massless fields, then picking up quadratic terms in the 
derivatives with respect to the vortex world-volume coordinates 

L = J dx'dx 2 C [Hrt(<Pi(ty)),W&{<pi(t,a?))] x=1 , (4.7) 

where <fi represents the set of moduli parameters (77, fj, 0, 0) or (a/, ft', a', V) contained in the 
moduli matrix f)3.20p . 
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If the coupling constant A is close to the BPS limit A = 1, we can still use the moduli approx- 
imation, to investigate dynamics of the non-BPS non-Abelian vortices by adding a potential of 
order |1 — A 2 | 1 to the massless sigma model 

L = J dx l dx 2 C [# sol (^(t,x 3 )), WU^(t,x 3 ))] x=1 - V{ifi). (4.8) 

We shall now calculate the effective potential V((pi) using the method suggested by Hindmarsh, 
who calculated this effective potential for non-BPS semilocal vortex in the Abelian-Higgs model 



First we write the Lagrangian ( 13. 2\\ in the following way 

I = Abps + (In ~ Hrff . (4.9) 

We get non-BPS corrections of order 0(A 2 — 1) by putting BPS solutions into Eq. ( 14.91) . This is 
because the first term is minimized by the BPS solution, while the second one is already a term 
of order 0(A 2 — 1). The energy functional thus takes the following form 

8 = 2^ = 2 + ^ 8^ ^ / dxldX " Tr i 1 ~ ^bps(^)^ps(^)) 2 (4-10) 

where i?BPs(v 9 «) stands for the BPS solution generated by the moduli matrices in Eq. (13. 20[) . The 
first term corresponds to the mass of two BPS vortices and the second term is the deviation from 
the BPS solutions which is nothing but the effective potential we want. 

In this section we consider the effective potential on the moduli space of coincident vortices. 
To this end, it suffices to consider only the following matrices 

Hg* = f ' ) , H™ = f \ ) . (4.11, 
\ z J \ —a z 1 I 

The parameters r\ and a! are related by r\ — 1/a'. The effective potential on the moduli space for 
two coincident vortices is thus I 2 ) 

V(\v\) _ (A 2 - 



2tTV 2 87T 



J dx l dx 2 Tr (l - #bps(M)#Jps(M)) 2 = (A 2 - l)V(hl), (4.12) 



where we have defined a reduced effective potential V which is independent of A. To evaluate this 
effective potential, we need to solve the BPS equations for a composite state of two non-Abelian 
vortices with an intermediate value of rj. Such numerical solutions found in [11]. We propose 
here another reliable technique for the numerics, which needs much less algebraic efforts. In the 



2 The potential depends only on \rj\ because the phase of rj can be absorbed by the flavor symmetry. The same 
holds for the coordinate a' 
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Magnetic Flux 




Fig. 3: Left: Numerical plots of Y\ (red), Y2 (green) and I3 (blue) for r/ = 3. The black line is Y\ = Y3 
for 77 = 0. Right: Magnetic flux Tri 7 ^ with 77 = (red), r/ = 3 (green) and 77 = 00 (a 1 = 0) (blue). 



moduli matrix formalism, what we should solve is only the master equation ( 13. 161) for Q with 
A = 1 and H(z,z) = H (z). Because of the axial symmetry of the composite vortex and the 
boundary condition at infinity: 

Q -h. H Q (z)Hl(z), (4.13) 

we can make a simple ansatz for Q. For example in the patch (1,1) we can write 

n( i,i) = ( Mr) -ve- i6 w 2 (r) \ ^ 
\ — 7]e %6 W2{r) Wz{r) J 

The advantage of the moduli matrix formalism is that only three functions Wi(r) are needed and 
the formalism itself is gauge invariant. Plugging the ansatz into Eq. (13.161) . after some algebra 
we get the following differential equations 

Y " + ;>? + = 1 - e " Yl(r2 + ^ (415) 

Y" + l Y ' - ni + rY{-rYi){l + rYi-rYj) _ _ 

1 UI2 H i „v/ „V\2 

-Y 3r 2 



1 - 


?7 2 e 5/1 


-Y 3 +2Y 2 


r 1:1 


|r7| 2 (H- 


rYi-rYif 


,y«2 | | 2 


- e Y!+Y 3 -2Y 2 



i-e'^r 2 , (4.17) 

where we have redefined the fields as u>j(r) = e ¥i ^ with i = 1,2,3. We solve numerically these 
differential equations using a simple relaxation method, see FigJSJ 

The effective potential can be obtained by plugging numerical solutions into Eq. (I4.12p . The 
result is shown in Fig. HI When we consider configurations with big values of \rj\, it is better 
to switch to the other patch and use the variable a' . We can still make a similar ansatz for Q, 
leading to simple differential equations like those in Eqs. fl4.15p ~ fl4.17p . In the left of Fig. S]we 
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0.01 0.02 0.03 0.04 



Fig. 4: Numerical plots of the effective reduced potential V(|r/|). In the second plot a' = 1/rj. 



show a numerical plot of the reduced effective potential V from \rj\ = (\a'\ = oo) to \r]\ = 20. 
In the right of Fig. H]we show another plot from \a'\ = (|^| = oo) to \a'\ = 1/20. To give an 
estimate of the range of validity of our approximation we can compare the results obtained in this 
section with the numerical integrations obtained for (2, 0)-vortices and also (1, l)-vortices. The 
comparison is made in Table [2]from which we can argue that the effective potential approximation 
gives result with an accuracy around 10% for the range of values 0.7 < A < 1.15. 



A 


(2, 0) num 


(2,0) eff 


(1; l)num 


(M)cff 


0.6 


1.52625 


1.65280 


1.62611 


1.73440 


0.7 


1.65337 


1.72332 


1.72880 


1.78835 


0.9 


1.88936 


1.89692 


1.91473 


1.92115 


0.95 


1.94523 


1.94711 


1.95793 


1.95954 


1 


2.00000 


2.00000 


2.00000 


2.00000 


1.05 


2.05376 


2.05561 


2.04102 


2.04254 


1.1 


2.10655 


2.11393 


2.08106 


2.08715 


1.15 


2.15843 


2.17496 


2.12018 


2.13384 


1.2 


2.20944 


2.23870 


2.15843 


2.18260 



Table 2: Numerical value for the masses of coincident vortices. (2,0) num is for the numerical results 
while (2, 0) e ff is for our approximation using the effective potential. 



In the type II case ((A 2 — 1) > 0) the effective potential has the same qualitative behavior as 
showed in the figure. As we expected, it has a minimum at \rj\ = 0. This matches the previous 
result that the (1, l)-vortex is energetically preferred to the (2, 0)-vortex. In the type I case 
((A 2 — 1) < 0) the shape of the effective potential can be obtained just by flipping the overall sign 
of the effective potential of the type II case. Then the effective potential always takes a negative 
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value, which is consistent with the fact that the masses of the type I vortices are less than that 
of the BPS vortices, see the middle column in Table [TJ Contrary to the type II case, the type 
I potential has a minimum at \a'\ = (\rj\ = oo). This means that the (2, 0)-vortex is preferred 
with respect to the (1,1) vortex. 

4.3 Interaction at generic vortex separation 

In this subsection we go on investigating the interactions of non-Abelian vortices in the U(2) 
gauge group at generic distances. As in Sect. 14.21 we will use the moduli space approximation, 
considering only small deviations from the BPS case. The generic configurations are described 
by the moduli matrices in Eq. (13.201) . We will consider here only the reduced (1, 1) patch defined 
in Eq. (13.211) . By putting the two vortices on the real axis we can reduce z to a real parameter 
d. Furthermore, by the flavor symmetry, we can freely put fj — and suppress the phase of rj. 
The relevant configurations will be described by the following moduli matrix: 




-10 -5 5 10 -10 -5 5 10 



Fig. 5: The magnetic flux Tri ? i2 of a configuration of 2 BPS vortices. The figure shows slices including 
the centers of the vortex, (red, green, blue) correspond to r/ = (0, 4, oo) while (solid, small broken, wide 
broken) lines to d = (0, 1, 4). The left panel shows the fine-tuned model with e = g and the right shows 
the model with e = 2g. 

Now let us study the effective potential as function of rj and d. We first need the numerical 
solution to the BPS master equation for two vortices with any relative distance and orientation. 
Unlike the computation for the coincident vortices, we do not have an axial symmetry. We 
can no longer reduce the problem to one spatial dimension by making an appropriate ansatz. 
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Nevertheless the moduli matrix formalism is a powerful tool also for the numerical calculations. 
The master equation is a 2 by 2 hermitian matrix, so it includes four real 2nd order partial 
differential equations. Despite the great complexity of this system of coupled equations, the 
relaxation method is very effective to solve the problem. We show several numerical solutions in 
Fig. [51 As before, once we get the numerical solution to the BPS equations, the effective potential 
is obtained by plugging them into Eq. fl4.10p . 

The numerical plot of the reduced effective potential V is shown in Fig. [61 The effective 
potential for the type II case has the same shape, up to a small positive factor (A 2 — 1). The 
potential forms a hill whose top is at (d, \r]\) = (0, oo). It clearly shows that two vortices feel 
repulsive forces, in both the real and internal space, for every distance and relative orientation. 
The minima of the potential has a fiat direction along the <i-axis where the orientations are anti- 
parallel (rj = 0) and along the t] axis at infinite distance (d = oo). Therefore the anti-parallel 
vortices do not interact. 




Fig. 6: Numerical plot of the reduced effective potential V(rj, d). 

In the type I case (A < 1) the effective potential is upside-down of that of the type II case. 
There is unique minimum of the potential at (d, \r]\) = (0, oo). This means that attractive 
force works not only for the distance in real space but also among the internal orientations. 
Configurations with anti-parallel orientations do not interact, but these configurations represent 
unstable points of equilibrium. Type I vortices always stick together. 

5 Vortices with generic couplings 

In this section we will shift from the fine-tuned U(N) model ( 12. II) to the more general model 
defined in Eqs. (12. 6ft and ( 12. 7ft . This will lead to more complicated algebra, but we will also 
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clarify interactions with different qualitative behaviors. There are 5 parameters (g, e, X g , A e , v) in 
the model, but we can reduce their number with the following rescaling 



H^vH, Wp^evWp, x^^. {5.1) 

ev 



Then the Lagrangian (I2.6P is expressed as follows 
t = Tr 



2„2 ^ 2 fJ,Vj 



2^ T rX 2 -^Tr(X°T°-l w ) 2 , (5.2) 



where £ = C/e v and we have introduced the ratio of the two gauge couplings 

7 = 9 -. (5.3) 

e 

We have three effective couplings 7, A e , \ g and the rescaled masses of particles are 

M u{1) = 1, M SU (N) =7, M s = A e , M ad = 7 A 9 . (5.4) 

The Lagrangian can be thought of as the bosonic part of a supersymmetric model only when 
both the parameters A e and X g are unity. In the BPS case we can easily find the BPS equations 

VH = 0, F l2 = -Ljt, F® 2 T° = l - (X°T° - 1 N ) . (5.5) 

The last two equations can be rewritten in a compact way 

Fn = ~ (X - 1 N ) + ^^X. (5.6) 

The first BPS equation in Eq. f)5.5p can be solved using the moduli matrix in the usual way 

H = S~ 1 (z, z)H (z), W = -iS- x dS, (5.7) 

where S is a GL(N, C) matrix. We would like to stress that this solution does not depend on 7 
so that the moduli space of the BPS vortices is the same as that of the well investigated vortices 
in the equal gauge coupling theory g = ejf| The Eq. ( 15.61) can be rewritten in a gauge invariant 
fashion as 



d (QdQ^) = \ (Q Q^ - l N ) + t-1 (q q^ - Tr( ^' 1} l^ 

where Q = is same as before and f2 = HqHq. 

3 The moduli space is the same from the topological point of view, while the metric will be different. 



(5.* 
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Now we are ready to investigate interactions between two almost BPS vortices by using same 
strategy as we used in Sec. 4. An effective action of the moduli dynamics for appropriately small 
|1 — Ag | <C 1 is obtained by plugging BPS solutions into the action. Then we get 



27TV 2 



Tt X 2 + ^LA Tr ( X °T° - 1 N \ 



— / dx 2 Tr 

2tt 



7T 

7 2 



8ti 



{F 12 y + (Xi-i)(F? 2 T°f 



(5.9) 



where we have used the BPS equations in the second line. Let us define the Abelian and the 
non-Abelian potentials as 



i0 rr0\2 



VeM; T )= / dx 2 Tr(F? 2 T 



V g (r],d-j) = / dx 2 Tr(F 12 



The true potential is a linear combination of them 



V(r), d- 7 , A e , \ g ) = (Ag - l)V e (v, d; 7) + 



7 2 



(5.10) 



(5.11) 



Notice that V e „ is determined by the BPS solutions, so it does not depend on A 



e,g- 



5.1 Equal gauge coupling 7 = 1 revisited 

In Sec. 4 we have discussed the effective potential for two vortices for any separation and with 
any relative orientation in a model with 7 = 1 and A = X g = A e . The potential is shown in Fig. [6l 
The effective potential comes from two pieces: the Abelian V e and the non-Abelian V g potential 
given in Eq. (15. 101) . In Figs. [7J and [8] we show V e and V g taking several slices of Fig. [61 

Let us consider the case with X 2 — 1 > and X 2 — 1 > 0. In this case the effective potential 
will have the same qualitative behaviors like the reduced potentials in the Figs. [7J [8] and 
The figures shows how V e and V g behaves very differently. In particular, the Abelian potential is 
always repulsive, both in the real and internal spacej (see the red dots in Figs. Hand E]). The non- 
Abelian potential is on the contrary sensitive on the orientations. In particular, Fig. [7] shows that 
it is repulsive for parallel vortices while it is attractive for anti-parallel ones. Furthermore, the 
non-Abelian potential becomes almost flat (along the spatial coordinate d) with the orientation 
77 ~ 4, see the middle of Fig. [71 The blue dots in Fig. [8] reveal that the non-Abelian potential 
always gives attractive forces in the internal space. When the two scalar couplings are equal, 
A 2 = X 2 , the left picture in Fig. [7J clearly shows how the two potentials exactly cancel for anti- 
parallel vortices, recovering the result of the previous section. 



1 Attraction in the internal space here just means that orientations tend to become the same. 
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anti-parallel intermediate parallel 

Fig. 7: Effective potential vs. distance of vortices d. Red dots denote V e while blue denote V g . All 
plots are for 7 = 1. The left figure shows the two anti-parallel vortices r/ = 0, the middle shows rj = 4.1 
and the right shows the two parallel vortices rj = 00. 
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d=0 d=l d=A 

Fig. 8: Effective potential vs. 77. Red dots denote V e while blue denote V g . 



Of course, the true effective potential depends on A e and X g through the combination in 
Eq. (15.111) . This indicates the interaction between non-Abelian vortices is quite rich in comparison 
with that of the ANO vortices. Let us make a further example. The rightmost panel in Fig. [7] 
shows that V e and V g for parallel vortices (77 = 00) are identical^. Thus if we take \ 2 e + X 2 = 2, 
the interactions between parallel vortices vanishes while between anti-parallel vortices they do 
not canceled out. This is just opposite to what we found in the case Ag = X 2 . We summarize the 
possible behaviors of the Abelian and non-Abelian forces, when the gauge couplings are equal 
e = 9 (7 = 1)) i n Table El 

5.2 Different gauge coupling 7 7^ 1 

We now consider interactions between non-Abelian vortices with different gauge coupling e 7^ g 
(7 7^ 1). Several numerical solutions are given in the left panel of Fig. [51 In Figs. [TU1 and [TTJ we 
show two numerical examples for the reduced effective potentials V e , V g given in Eq. (15.101) . 

The plots show that the qualitative features of V e and V g are basically the same as what is 
5 This statement can be proved analytically. 
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d 



4 

d 



Fig. 9: The Abelian potential V e (left) and the non-Abelian potential V g (right) for 7 = 1. 





A 2 >1 


*J=1 




XI >1 




N=0 
A=+ 




A e 2 = l 




N=0 
A=0 




A e 2 <l 


N = (-, + ) 


N=0 
A=- 


N = ( + ,-) 



Table 3: The forces between two vortices in the 7 = 1 case. N and A stand for non-Abelian force and 
Abelian force, respectively. means no force, + means repulsive and — means attractive. N=(— ,+) 
means that there is an attractive force for anti-parallel vortices and a repulsive force for parallel vortices. 

discussed in the equal gauge coupling case (7 = 1). Therefore, the qualitative classification of 
the forces given in Table [3] is still valid for 7 7^ 1. We observe that the Abelian potential tends, 
at large distances, to a value smaller than the non-Abelian one for 7 < 1, while opposite happens 
for 7 > 1. The only things that have non dependence on 7 are the values of the potentials at 
(1, l)-vortices with (d,rj) = (0,0). Regardless of the gauge couplings V g (r) — 0, d — 0) = while 
V e (r] = 0, d = 0) '= 0.41. This is because the corresponding solution is proportional to the unit 
matrix, so that there are no contributions from the non-Abelian part. 

The effective potential is obtained from the linear combination in Eq. ( 15. lip and also depends 
on three parameters 7, A e and X g . With this big freedom we can obtain a lot of interesting 
interactions. For example, we can have potentials which develop a global minimum at some 
finite non zero distance. In such cases two vortices may be bounded at that distance. We can 
show a concrete potential in Fig. [121 The figure shows the presence of a minimum around d ~ 2^ 
This kind of behavior have not been found for the ANO type I /II vortices and the possibility of 

6 These kind of potentials are really generic. In fact, the presence of this minima docs not require strong 
constraints on the couplings. 
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Fig. 10: Effective potential with 7 = 1/2 vs. separation, (red, blue) = (V e , V g ). 
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Fig. 11: Effective potential with 7 = 1.3 vs. separation, (red, blue) = (V e , V g ). 



bounded vortices really results from the non-Abelian symmetry. 



6 Interaction at large vortex separation 
6.1 Vortices in fine-tuned models e = g and A e = X g 

In this subsection we will obtain an analytic formula for the asymptotic forces between vortices 
at large separation. We follow the technique developed in Refs. [HHl HSj . First of all, we need to 
find asymptotic behaviors of the scalar and the gauge fields. We again consider the (1, 0)-vortex 

*««-(;:). *--(;)■ («) 

When we are sufficiently far from the core of the vortex, Y(r) and f(r) in Eqs. (14. 2 j) and (14.31) 
can be written as 

Y = 2\ogr + 5Y, f = l + 5f (6.2) 
7 Similar behaviors in the static intervortex potential were also observed for Z3 vortices in [37| . 
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Fig. 12: 7 = 1/2, A e = 1.2, X g = 1.06: From rj = (green) to r) = 7 (blue) with d = ~ 5 for each 77. 



where <5F and Sf are small quantities. Plugging these into Eqs. (14. 2p and (14. 3 p and taking only 
linear terms in SY and Sf, we obtain the following linearized equations 

5Y" + -SY' - X 2 5Y - 2 f Sf" + -Sf - X 2 5f] = 0, (6.3) 
r \ r J 

SY'" + -<jy" - —SY' - 5Y' = 0. (6.4) 
Solutions to these equations are analytically obtained to be 

SY - 25 f = -±K {\r), SY = -K (r) - C, (6.5) 

7T 7T 

where A (r) is the modified Bessel's function of zeroth order and q, m and C are integration 
constants. C must be because SY — » as r — > 00 while q, m should be determined by the 
original equation of motion. In the BPS case, Sf = (/ = 1), thus q = —m. Eq. (I6.5P leads to 
the well known asymptotic behavior of the ANO vortex 

= fe^ Y z ~(l + 5f- V) e w = (l + ^o(Ar)) e ie , (6.6) 

%i] = 4aF^4e^(21ogr + ^) = -l(i-£AMr))e- (6.7) 

where A x = — Aq and we have defined #[1,1] and Wp.,1] as [1, 1] elements of H and W in Eq. (13.71) 
with the k = 1 ansatz (14.11) . 

Next we treat the vortices as point particles in a linear field theory coupled with a scalar 
source p and a vector current j M . To linearize the Yang-Mills-Higgs Lagrangian, we choose a 
gauge such that the Higgs fields is given by the following hermitian matrix 











: 





h° + h 3 h 1 — ih 2 \ if w° + wf, wj. - 



h=\:: i + ^i ■;„ i, w, = ^[ 7;;^ 7 -7 |. (6. 



/i 1 + ih 2 h° — h 3 / 2 \ y,i _)_ f w 2 _ ^3 
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Here all the fields h a , w" are real. Then the quadratic part of the Lagrangian ( 13. 2p is of the form 



(2) 
free 



E 

a=0 



1 



1 



A 2 



a\2 



4 



(6.9) 



where we have defined the Abelian field strength /" = d^wf, — d v w' \. We also take into account 
the external source terms to realize the point vortex 

3 

Source = ^2[ P a h a -j^]. (6.10) 

a=0 

The scalar and the vector sources should be determined so that the asymptotic behavior of the 
fields in Eqs. (16. 6p and (16.71) are replicated. Equations of motions are of the form 



{n + x 2 )h a = P a , (□ + !)< = j;. 



(6.11) 



In order to replicate the (1, 0)-vortex corresponding to the k = 1 ansatz (14.11) . we just need to 
mimic the result of Refs. [381 H3] because the single non- Abelian vortex is a mere embedding of 
the ANO vortex as mentioned earlier. In fact, only (h° , w®, p° , j®) = {h 3 , w 3 , p 3 , j 3 ) are relevant 
and all the others are zero: 



h° = h 3 



2vr 



K (Xr), p°=p 3 = q 8(r), 



w° = w 3 



m 
2^ 



^kxV/f„(r), 



j°=f 



-mk x V5(r) 



(6.12) 



where k is a spatial fictitious unit vector along the vortex world- volume. The vortex configuration 
with general orientation is also treated easily, since the origin of the orientation is the Nambu- 
Goldstone mode associated with the broken SU (2) color-flavor symmetry 

a* \ ( a 8 

8* J \ -8* a* 

where \a\ 2 + \8\ 2 = 1. The fields H and receive the following transformations, keeping the 
hermitian form of (16. 8p . 



(6.13) 



X 




ui 



X 




\a\ 2 a*8 
ad* \8\ 2 



X. 



(6.14) 



The scalar interaction between a vortex at x = x x with the orientation (f>i and another vortex at 
x = x 2 with the orientation 2 can be obtained by 

/i°(x-xi) 




dx 2 Tr 



\a\ 2 a*8 
ad* \8\ 2 



P (x- x 2 J 



a 



2tt 



K (A|x 1 -x 2 |). 



(6.15) 
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The gauge interaction is also obtained by similar way 

xi) 



dx 2 Tr 



w°fx 








a 



a*(3 



a(3* \(3\' 



,m 



-\a\ —K (\xx - x 2 |). 



2tt 



(6.16) 



Then total potential is V\ 



int 



¥2 



-L h 
2 



mi 





2 




$1 







— K (X\x.i - x 2 | 



^A' (|xi -x 2 | 



(6.17) 



where lal 2 



0102 



■2 -2 is invariant under the global color-flavor rotation. When two vortices have 

I 01 1 1 02 1 

parallel orientations, this potential becomes that of two ANO vortices [3H]. On the other hand, 
the potential vanishes when their orientations are anti-parallel. This agrees with the numerical 
result found in the previous sections. In the BPS limit A = 1 (q = m), the interaction becomes 
precisely zero. 



Since A (Ar) ~ y/ir /2Xre Ar , the potential asymptotically reduces to 

~ Xr for A < 1 Type I 
for A > 1 Type II 





|0l0 


I 2 a 2 




01 


•It 
- 1 


2 |"27r 
m 2 


|0l 


|2 


(f>2 


\ 2 2n\ 



IX 

2A7' 



7T 




(6.18) 



2r 



where r = |xi — x 2 | ^> 1. If we fix the relative orientation being some finite value, the force 
F r = — d r Vi nt between two vortices is attractive for A < 1 and repulsive for A > 1 similar to 
the force between ANO vortices. The force vanishes when the relative orientation becomes anti- 
parallel. If we fix the distance by hand, the orientations tend to be anti-parallel for the type II 
while the parallel configuration are preferred for the type I case. 



6.2 Vortices with general couplings 

It is quite straightforward to generalize the results of the previous section to the case of generic 
couplings. We can of course use the same gauge as in Eq. (]6.8p . The quadratic Lagrangian (12.61) 
is of the form 
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The external sources can be still reproduced by source terms as in Eqs. (16.101) . ( 16.111) . The 
linearized equations following from the above Lagrangian are of the form 

(^□ + i)< = #, (n + xy)h a = P a , (p + i)wl = j% (n + \l)h° = p°. (6.20) 

For the (1, 0)-vortex the only non-zero profile functions are (h°, w®, p°, j®) and (h 3 , w 3 , p 3 , j 3 ). 
But these profiles are no longer equal and we need to deal with them independently. The only 
difference from the similar equations (16.111) is for the masses of the particles. The masses are 
directly related to asymptotic tails of vector and scalar fields. We can easily find the solutions 
by doubling Eqs. ( 16.121) 



h = ^K (X e r), w° = -^kxVK (r), p° = q °S(r), j° = -m° k x V5(r), (6.21] 
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h? = 1-K {\ gl r), w 3 = -— k x VK (jr), p 3 = q 3 S(r), j 3 

Z7l ZTX 



-m 3 kxV5(r). (6.22) 



The vortex with the orientation <p2 in Eq. (16. 131) can be obtained by performing an 5'[/(2)c+f 
rotation like Eq. (I6.14p 
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Similar to Eqs. (16.151) and (I6.16p . we find the total potential Vi n t 
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When we tune the parameters to 7 = 1 and X g = X e (q° = q 3 , m° = m 3 ), this effective potential 
is exactly identical to that of Eq. ( 16.171) . In the BPS limit A e = X g — 1, the interaction becomes 
precisely zero because q° = m° and q 3 = m 3 . 

At large distance, the interactions between vortices are dominated by the particles with the 
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lowest mass M\ ow . There are four possible regimes 
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which exhaust all the possible kinds of asymptotic potentials of this system. This generalizes the 
type I/II classification of Abelian superconductors. We have found two new categories, called 
type I* and type II*, in which the force can be attractive or repulsive depending on the relative 
orientation. The type I force is always attractive and the type II force is repulsive regardless of the 
relative orientation. On the other hand, type I* and type II* depend on the relative orientation. 
In the type I* case the forces between parallel vortices are attractive while anti-parallel vortices 
repel each other. The type II* vortices feel opposite forces to the type l1j. Note that we have 
used the same terms type I/II for the fine-tuned model in sect. I6.ll In the perspective of this 
section, they should be called type I+F/II+II* because of the degeneracy of some masses. It is 
interesting to compare these results with the recently studied asymptotic interactions between 
non-BPS non- Abelian global vortices [23]. A very different feature of global vortices is that the 
interactions are always repulsive. This is because they are mediated by Nambu-Goldstone zero 
modes, whereas in our model these particles are all eaten by the gauge bosons thanks to the 
Higgs mechanism. 

We find a nice matching of qualitative features between the numerical results of the previous 
sections and the semi-analytical results in this section. Let us look at Figs. [10] and [TTJ In 
all the cases, we found that the Abelian potentials are attractive regardless of the orientations 
while the non-Abelian potentials are sensitive to those. These properties are well shown also in 
the semi-analytical results in f|6.25p . The type I/II interactions originated by the U(l) part are 
independent of the orientations whereas the type I*/II* which are coming from the SU(N) part 
do depend on them. 

The result in Eq. ( 16.251) is easily extended to the general case of U(l) x SU(N). This can be 
done by just thinking of the orientation vectors as taking values in CP^" 1 . 



8 In Ref. [41j we have found a similar result in a supersymmetric theory. The relation between the two notations 
is explained in the Appendix. The supersymmetric theory in Ref. [H] shows only type I/I* behaviors. 
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7 Conclusion and discussion 



In this paper we have studied static interactions between non-BPS vortices in SU(N) x U(l) 
gauge theories with Higgs fields in the fundamental representation. We have discussed models 
with arbitrary gauge and scalar couplings. We have numerically computed the effective potential 
for almost BPS configurations for arbitrary separations and any internal orientations. We have 
also obtained analytic expressions for the static forces between well separated non-BPS vortices. 
This expression is valid also for models far from BPS limit. 

For the fine-tuned model we found interaction pattern similar to that of the ANO vortices 
in the Abelian-Higgs model. The numerical effective potential is given in Fig. El it depends on 
both the relative distance and the orientations. The asymptotic potential between two vortices is 
given by Eq. (16.18p . In this model the mass of the U (1) and of the SU(N) vector bosons are same 
Msu(N) — Mum, and also all the scalars have the same masses M s = M^. We thus have only 
two mass scales, which corresponds to two different asymptotic regimes. For A < 1 (M s < Mu(i)) 
there is universal attraction (type I) and for A > 1 (M s > Mum) universal repulsion (type II). 
Both the numerical and the analytical result show that the interactions between two anti-parallel 
vortices vanish; this configuration is unstable for type I vortices and stable for type II. So in 
this last case the part of the moduli space which corresponds to vortices with opposite CP 1 
orientations at arbitrary distance survives the non-BPS perturbation. 

In models with arbitrary couplings, on the other hand, the pattern of interactions becomes 
richer. In this case we considered separately the Abelian and non-Abelian contributions V e 
and V g to the effective potential. The two show very different qualitative behavior. While the 
Abelian contribution is always attractive (or repulsive) for a given choice of the parameters, the 
non-Abelian one can be attractive or repulsive depending on the relative internal orientation of 
the two vortices. These properties combined with the fact that the full effective potential is the 
linear combination given in Eq. (15. lip deduce that we can obtain many qualitatively different 
type of interactions depending on the choice of the parameters of the theory. Such a variety 
also appears in the possible asymptotic behavior of the interactions. In the theory there are four 
different mass scales, the masses of the SU(N) vector bosons Msu(n)> the U(l) vector boson 
M[/(i), the adjoint scalars M a d and the singlet scalar M s under the color-flavor symmetry. This 
leads to four different asymptotic regimes, classified in Eq. (16.251) . We found, in addition to type 
I and type II, new types of interaction mediated by the non-Abelian particles, which we call type 
I* and type II*. The type I (type II) force is attractive (repulsive), and occurs when the singlet 
scalar (U(l) vector boson) has the smallest mass. These forces do not depend on the relative 
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orientation. In the type I* case there is an attractive force for parallel orientations and repulsive 
for anti-parallel ones where the asymptotic force is mediated by the lightest non-Abelian scalar 
fields. On the other hand for type II* there is repulsion for parallel orientation and attraction for 
anti-parallel ones and the force at large distance is mediated by the lightest non-Abelian vector 
field. 

The dynamics of the interactions of the non-BPS non-Abelian vortices is quite rich. Let us 
give comments on some possible further directions: 

• Reconnection rate of the cosmic string: The slow moving non-Abelian BPS vortex strings, 
as cosmic strings, were shown to always reconnect with probability one [T7J [18]. This is 
important in order to distinguish solitonic cosmic strings from fundamental cosmic strings 
which generically have very small reconnection probability. Consider the two dimensional 
space spanned by 77 and d; the coordinate d corresponds to the vortices relative distance and 
rj to the internal colour-flavour orientation. The (2, 0)-vortices are at rj — > 00 and the (1, 1)- 
vortices at rj = 0. Only the orbits which pass through the point d = 0, 77 = correspond 
to scatterings where no reconnection occurs; these orbits represent scatterings with very 
finely tuned initial conditions. The fine-tuned scattering process cannot contribute to the 
reconnection probability; for this reason the reconnection rate of the non-Abelian BPS 
vortex is one [T71 [18] . 

When we consider strong non-BPS corrections, the moduli space approximation is no longer 
valid and this conclusion could drastically change. For some values of the couplings there will 
appear regimes in which the coincident (1, l)-vortices are favored to the coincident (2,0)- 
vortices. In that case we may expect that the reconnection probability becomes smaller 
than one. As a very simple example, we can consider the fine tuned model for A > 1. In 
this case the energetically favored configuration is the one with two vortices with opposite 
CP 1 orientation, which never reconnect. In this case we can expect a reduction of the 
reconnection rate. On the contrary for A < 1, when the (2, 0) state is energetically favored, 
we expect that the reconnection rate should still be one. It would be interesting to make a 
detailed numerical investigation of the scattering process of non-BPS non-Abelian vortices, 
in order to clarify how the non-BPS corrections could modify the reconnection rate. 

• Non-Abelian vortices and Abrikosov lattice: In the usual type II superconductor (the Abelian- 
Higgs model), if a large number of vortices penetrate a region of given area A, they will form 
a hexagonal lattice (Abrikosov lattice) rather than forming a square one [39] . This is verified 
experimentally. If we consider the same for non-Abelian vortices, we expect that the prop- 
erty of the lattice can be quite different, say, lattice spacing and/or form can change. An 
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interesting possibility is the appearance of phase transitions due to the change of the lattice 
structure when the density of vortices change. This eventuality is suggested by the possible 
presence of interactions which change with distances. For example in the case of Fig. [T2l 
there are repulsive forces at short distances while at large distances they are attractive. 

• Quantum aspects: In this paper we focused on the classical aspects of the interactions 
between non-Abelian vortices. In the theoretical set-up that we have discussed, the quantum 
aspects of the infrared physics of a single vortex are described by an effective bosonic CP N ~ l 
sigma model (the theoretical setting is in this sense similar to the one discussed Ref. [7]). 
Let us consider two vortices at large distance; in the type I* and in the type II* regimes 
the quantum physics will be described by two CP N ~ X sigma models with an interaction 
potential given by Eq. (16.251) . It would be interesting to study the effect of this term in 
the sigma model physics. Another interesting problem is the numerical determination of 
the effective theory for vortices at generic separations. To do this one has to determine the 
Manton metric on the full moduli space. 
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A Note on the relation between two formalisms 

The aim of this small section is to relate the moduli matrix formalism that we used in this 
paper to the direct ansatz approach we choose in a previous work [JT]. To this end we must 
find the relation between the angle a used in jH] and the moduli matrix parameters rj = 1/a'. 
The resulting relation between the two variables, in the case of coincident vortices, is rather non 
trivial. 

Coincident vortices. In the moduli matrix formalism, two coincident vortices are described by 
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the following moduli matrix 

*r =(;-:), 

while in [JT] we used explicit ansatzs for the fields. In particular, for the squarks field we 
have used 

/ -cosfe 2 ^«i(r) sinfe^/caM \ 

#ans = • (A.2) 

\ — sin |e l</3 K 3 (r) — cos^K^r) J 

This form for the squark fields lead us, in [HI [12], to conjecture that this solution is asso- 
ciated with that given by the following moduli matrix: 

(— cos |z 2 sin | z \ 
(A.3) 
— sin |z — cos | J 

Using a V equivalence we can put H ans on the standard (1, 1) form 

nr = vwm^ = ( * cot( 2 a/2) ) . (A.4) 

This simple argument gave us the relation: rj = cot (a/2). 

Here we point out that this conjecture is wrong. In fact the functions K 2 , ^3, ^4) in 
Eq. (1A.2|) have an implicit dependence on a, which come out only after solving the differential 
equations for the vortices. This means that it is difficult to find a simple analytical relation 
between the parameters of the two formalisms. Also, this relation is not the same for all the 
values of 7 = g/e. In Fig. [13] we plotted the numerical relation between a and r] for 7=1. 
The corresponding values for the two variables are found making an empirical match of some 
gauge invariant functions which have a non-trivial dependence on the relative orientation. 
One of these functions, TrF 12 , is showed in Fig. [31 The fact that a perfect matching is 
possible between gauge invariant functions in both the formalisms is a strong numerical 
check of the consistence of the two approaches. 

Well separated vortices. Vortices at large separation have well-defined global orientations in 
the internal space, and thus a well-defined notion of relative orientation. This enables us 
to find an exact relation between the parameters of the two formalisms. In [?T] the global 
orientation of a non-Abelian vortex was defined by a vector which takes values in the internal 
space: n. Using global color-flavor transformations we can always put the orientations n\ 
and n 2 of two vortices on the following standard form 

fii = (0,0,1), n 2 = (— sin a, 0, cos a). (A. 5) 
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Fig. 13: The solid line is the naive relation 77 = cot (a/2) and the dots are the true numerical relations 
between 77 and a for 7 = 1 



The vector n 2 can be obtained from fii acting with a global rotation 
n 2 ■ f = U^r 3 U, with U = 



cos y — sin g 



sin I cos I 



(A.6) 



We can repeat the same argument for the orientation vectors defined within the moduli 
matrix formalism (see Eq. (13.181) ). We can choose the following identification: 



(fii, n 2 ) = ((0, 0, 1), (— sina, 0, cosa)) <^ 
If we act with the same global rotation on 0i 

02 = rffl 
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we find the relation b' = — tan f . With this identification we can easily check the identity: 
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(A.9) 



This shows the consistence of the expressions for the asymptotic forces obtained in this 
paper in Eq. (I6.25P with those found in 
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